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ABSTRACT 

We present ALMA continuum and spectral line observations of the young Brown Dwarf p-Oph 102 at 
about 0.89 mm and 3.2 mm. We detect dust emission from the disk at these wavelengths and derive 
an upper limit on the radius of the dusty disk of ~ 40 AU. The derived variation of the dust opacity 
with frequency in the mm provides evidence for the presence of mm-sized grains in the disk outer 
regions. This result demonstrates that mm- grains are found even in the low density environments of 
Brown Dwarf disks and challenges our current understanding of dust evolution in disks. The CO map 
at 345 GHz clearly reveals molecular gas emission at the location of the Brown Dwarf, indicating a 
gas-rich disk as typically found for disks surrounding young pre-Main Sequence stars. We derive a 
disk mass of ~ 0.3 — 1% of the mass of the central Brown Dwarf, similar to the typical values found 
for disks around more massive young stars. 

Subject headings: circumstellar matter — brown dwarfs — stars: individual (p— Oph 102) — planets 
and satellites: formation — submillimeter: stars 



1. INTRODUCTION 



Brown Dwarfs (BDs) are very low-mass 
(M<O.O85M ), stellar-like objects unable to burn 
hydrogen. It is now established that y oung BDs are 
surrounded by dusty circumstellar disks (|Natta fc Testi 
" INatta et al.| [2002] |Testi et aL] |20(^| l^choIzT 
ardhana 2008 ) and undergo a T Tauri-like phase 



Jayawari 

during their e arly evolution that involves disk accretion 



and outflows (I Jayawardhana et al. 2003 Natta et al. 



2004]). 



"Because of their very low masses, disks around BDs 
also represent an extreme environment to test planet for- 
mation theories. There are a few wide orbit, giant plan- 
ets known as companions to BDs, most likely formed 
as a result of the fragmentation of the proto-BD core 
or by gravitational instabilities in the outer disk (e.g 
2M1207b, 2M0441b, ~" " " " 



Chauvin et al.|[2005 Lodato et al. 



2005 Todorov et al. 



2010^ However, most extrasolar 



planets (and all rocky pla nets) are thought to form via 
the core accretion process ( |Matsuo et al | |2007), and ev- 
idence of the first steps of this process has been found 
in disks around more massive pre-Main Sequence (PMS) 
stars. In particular, sub-mm observations, which probe 
the bulk of the disk cold material, have shown that the 
grains h ave grown to to reach millimeter and centimeter 
size (e.g. Beckwith & Sargent 1991; jRodmann et al.|2 006; 
|Kicci et al ||2010al|Ubach et al.||201^| ~ 

In BD disks, infrared spectroscopy has shown that in 
several sources silicates on the disk surface have been 
significantly modified in size and crystallinity, as in PMS 
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stars ( jSterzik et aLl|2004| |Apai et al.||2005| ). Until now, 
due to the limited sensitivity of sub-mm telescopes, the 
detection of long wavelength emission from young BDs 
has been limited to co ntinuum photometr y of the few 



'Y 

brightest objects (e.g. |Scholz et al.||2006| [Mohanty et 



|al.||2012| . The start of ALMA operations allows us to 
perform more detailed studies. Here we report on multi- 
wavelength continuum and CO( J = 3 — 2) observations 
on the young BD p-Oph 102. This object (~ 60 M Jnp , 
M6-spectral type) in the p-Ophiuchi star forming region 
(SFR) is known to be surrounded by a disk from in- 
frared observations, to have a significant mass acc retion 
rate and to drive a wind and molecular outflow ( Bon- 



temps et aL [200lUNatta et al. 



2004; Whelan et at 



2005tPhan -Bao et al.||2008| ). In this Letter we report on 
our ALMA observations, which were designed to derive 
a solid estimate of the mm spectral index to constrain 
the dust properties. In particular we found evidence for 
grain growth to mm-grain sizes in the outer regions of 
the p-Oph 102 disk. Also, we clearly detected molecular 
gas emission from the disk and derived an estimate for 
its mass. 

2. ALMA OBSERVATIONS AND RESULTS 

2.1. Observations and data reduction 

We observed p-Oph 102 using ALMA Early Science in 
Cycle at Band 7 and 3 (about 345 and 100 GHz, respec- 
tively). Observations in Band 7 were performed using 15 
antennas in the Compact array configuration on Novem- 
ber 3, 2011, and 16 antennas in the Extended array con- 
figuration on May 23, 2012 (projected baseline lengths 
in the range from the shadowing limit to ~400 m). Ob- 
servations in Band 3 were conducted with 16 antennas 
in the Extended configuration on May 8, 2012 (projected 
baseline lengths from ~35 to ~402 m). All observations 
were done in good and stable weather conditions, with 
precipitable water vapor of ~ 0.8 — 1.2 mm and ~ 1.9 mm 
at Band 7 and 3, respectively. The ALMA correlator was 
set to record dual polarization with four separate spec- 
tral windows, each providing a bandwidth of 1.875 GHz 
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with channels of 0.488 MHz width. Spectral windows 
were centered at 331.103, 332.998, 343.103, 344.998 GHz 
for Band 7 observations and 86.1021, 87.9976, 98.103, 
99.999 GHz for Band 3. The total integration times 
on p-Oph 102 were approximately 15 and 30 minutes 
in Band 7 and 3, respectively. 

The interferometri c visibility data were r educed using 
the CAS A package QMcMullin et al.||2007[ ). J1924-292 
and 3C279 were observed as bandpass calibrators, Nep- 
tune and Titan for flux calibration. Simultaneous obser- 
vations of the 183 GHz water line with the water vapour 
radiometers were used to reduce atmospheric phase noise 
before using J1625-254 for standard complex gains cal- 
ibration. The flux scale was tied to the Butler- JPL- 
Horizons 2010 models of Neptune and Titan, resulting 
in an accuracy of ~ 10%. 

Imaging of the calibrated visibilities was done in 
CASA. For each band, a continuum map was produced 
using natural weighting and combining all the channels 
without line emission. The only line emission detected 
was the CO ( J = 3 — 2) rotational line in Band 7 which 
was imaged separately. The CO (3-2) shows significant 
extended emission as well as compact structure. We im- 
aged the CO data in two separate ways to highlight the 
different components: we first used all the data, natural 
weighting and a gaussian taper of the visibilities to re- 
cover as much as possible of the extended emissi on and 
compare with th e low angular resolution map of |Phan-] 



ity amplitude with projected baseline length, neither at 
the shorter nor at the longer baseline lengths probed by 
our observations. We conclude that our observations did 
not spatially resolve the disk continuum emission. Con- 
sidering our angular resolution of ~ 0.6" at 0.89 mm, 



Bao et al. ( |2008[ ); additionally, we also imaged the data 
from baselines exceeding ^ 60 kA with robust weighting 
to filter out emission on scales larger than about 3", or 
^ 400 AU at the Ophiuchus distance, and highlight the 
compact emission associated with p-Oph 102. 

2.2. Continuum maps 

Figure [I] shows the continuum maps at about 0.89 mm 
(Band 7) and 3.2 mm (Band 3). Dust emission from the 
disk surrounding p-Oph 102 is clearly detected in both 
maps. We measured flux densities of 4.10 ± 0.22 mJy at 
0.89 mm and of 0.22±0.03 mJy at 3.2 mm. The obtained 
rms-noise levels are consistent with the theoretical expec- 
tations given the on-source integration times. This gives 
a sub- mm spectral index a (F u oc v a ) between 0.89 and 
3.2 mm of «o. 89-3. 2mm = 2.29 ± 0.16, where the uncer- 
tainty accounts for both the flux calibration and the rms 
noise in the maps. This value lies well within the range 
of values measured for disks around single PMS stars in 
the p-Oph and Taurus SFRs (see Figure|3|. This suggests 
that the physical properties of the dust grains emitting 
in the sub-mm, e.g. their size, do not significantly vary 
in disks around PMS stars and BDs. The implications 
of this result will be discussed in Section 03 

Usin g observations with the Submillime ter array 



(SMA, [Ho et al.|[2004| |Phan-Bao et al | (l2008| ) reported 
a low signal-to-noise flux density ol 7 ± 3 mJy at about 
1.3 mm. Interpolating between our two measurements, 
we expect a flux density of 1.7 ± 0.3 mJy at 1.3 mm, 
lower than the value previously reported, although the 
discrepancy is only at the 2a-level. Note that already 
with ~ 15 minutes on-source with ALMA in Early Sci- 
ence we reached a sensitivity ~ 10 times better than the 
SMA. 

By inspecting the visibility data, within the uncertain- 
ties, we did not find evidence for a decrease of the visibil- 



an d a di s tance of 130 pc for the p-Oph SFR (Wilking et 
[aT]|2008| |Lombardi et al.|[2008[ ), this indicates that the 
dust is concentrated within less than ~ 40 AU from the 
central BD. 

2.3. CO molecular gas map 

The low angular resolution, tapered natural weight- 
ing map shows a morphology of the CO (J = 3 — 2) 
that closely resembles the map of the CO (J = 2 — 1) 
emission presented by |Phan-Bao et ah ( 2008| ). We con- 
firm the same main emission features as observed in the 
SMA map. They interpreted these structures as a bipo- 
lar molecular outflow emitted by p Oph 102. Even if the 
signal to noise in our map is significantly better than in 
the SMA data, our data are most likely significantly af- 
fected by missing flux at short baselines. Furthermore, 
the most prominent features attributed to the outflow are 
close to the edge of the ALMA primary beam, making it 
difficult to accurately determine the emission and its ex- 
tent. For these reasons, we will not discuss the extended 
emission and outflow further in this Letter. 

In the CO (J = 3 — 2) low angular resolution map we 
also detected compact emission coin ciden t with the con- 
tinuum emission described in Sect. |2.2 Left panel of 



Fig. [2] shows the map of integrated CO( J = 3 — 2) emis- 
sion, obtained by filtering out the extended CO emis- 
sion as described at the end of Section [2jl. Compact 
emission, that we associate with p-Oph 102, is clearly 
detected in several channels. The disk is detected at 
a > 10a level, with a velocity-integrated total flux of 
about 530 ± 45 mJy-km/s. Right panel of Fig. [2fehows 
the intensity weighted velocity field of CO molecular gas. 
We tentatively find the presence of a velocity gradient 
along the disk spanning a range of radial velocities of 
& 1 km/s. This is consistent with gas in Keplerian ro- 
tation in an inclined disk orbiting at distances > 10 AU 
from a 0.06 M BD. However, the angular (~ 0.5") and 
velocity (~ 0.4 km/s) resolutions of our observations do 
not allow a full characterization of the rotation curve of 
the disk and non-Keplerian rotation curves are still pos- 
sible. 

3. DISCUSSION 

To derive the disk mass and physical properties of the 
grains we model the mm continuum emission using a two- 
layer (surface+midplane) model of a flared disk heated 
by the radiation of the central BD (Chiang & Goldreich 
1997| |Dullemond et al.|2001[) . 

The adopted sub- stellar photospheric properties of p- 
Oph 102 and the dust model are described in the caption 
of Figure [4j Since our observations did not allow us to 
spatially resolve the dust emission of the disk, we could 
not strongly constrain the disk radius and inclination. 
Our analysis of the unres olved disk photo metry follows 
the procedure outlined in Testi et al. (2001), adapted for 
the case of young BDs. For the disk surface density we 
adopted a radial profile S oc r" p with p = 1, which is the 
typical value found for T Tauri disks in the sub-mm (e.g. 
|Isella et al.1|2009l |Andrews et~aT1[20Tol |Guilloteau etaT] 
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20TT|. A large value of the disk in clination is suggested 



by the detection of an optical jet (Whelan et al. 2005) 
and by the nearly symme tric morphology of the bipolar 
outflow in molecular CO ( |Phan-Bao et al.||2008 ). At the 
same time, the visual extinction of the~ 7^Uph 102 pho- 
tosphere, i.e. Ay ~ 3 ( |Natta et al.||2002[ ) poses an upper 
limit of ~ 80 degrees for the disk inclination, as a more in- 
clined disk would absorb t he (sub)-stellar radiation much 
more efficiently (see Skemer et al.||20lT ). We assume in 
the following an inclination of 70 degrees; however, our 
results do not depend critically on this choice. 

3.1. Disk mass and size 

The main results of our modeling can be derived from 
Figure [4j which shows the ALMA data overlaid with the 
predictions from the disk models presented above. Disk 
models with radii larger than > 5 AU can reproduce the 
data only if the spectral index of the dust opacity j3 {k v oc 
is smaller than 1. This is because for these disks the 
dust emission is predominantly coming from the optically 
thin outer regions of the disk. In this case, the spectral 
index of the SED reflects the dust opacity spectral index, 
and a relatively low value of f3 is needed to explain the 
measured spectral index ao. 89-3. 2mm- Models with disk 
radii < 5 AU reproduce the observed continuum flux and 
spectral index with a broad range of j3— values. This is 
because very small disks have to be relatively dense to 
explain the observed sub-mm fluxes: the dust emission 
in these models becomes mostly optically thick and, in 
this regime, it does n ot depend on the dust opacity law. 
(cf. the discussion in |Testi et al.||200l| ). 

Using the CO (J = 3 — 2) detection we can test the 
hypothesis of a small, optically thick disk. The mean 
brightness temperature of the CO (J = 3 — 2) in our 
- 0.5" FWHM beam is -8.3 K. As the gas in the disk 
is primarily heated by radiation from the central BD, 
its temperature cannot be very high, i.e. < 30-40 K, 
which in turns means that the beam filling factor of the 
disk cannot be very small. Under the assumptions of 
optically thick gas emission and face-on geometry, this 
argument sets a lower-limit of about 15 AU for the disk 
radius. For example, if the disk radius was 5 AU the 
beam filling factor would be < 2.3%, the exact value de- 
pending on the disk inclination. This would correspond 
to a gas temperature > 400 K, which is much higher than 
expected for a BD disk. 

The /3-values constrained by our analysis for these 
larger disks, i.e. 3 « 0.4 — 0.6, correspond to dust opac- 
ities values of % 89mm = ^ — 4.5 cm 2 /g (see also 
Ricci et al. 2010a). Given the measured flux densities, 



we derive estimates for the dust mass of ~ 2— 6-10 -6 M , 
or — 2 — 6 • 10 -4 M assuming an ISM-like gas-to-dust 
mass ratio of 100. These values correspond to a disk 
which contains —0.3-1% of the mass of the BD+disk sys- 
tem and are consistent with the distribution of values fo r 
more massive T Tauri systems ( Williams fc Cieza|2011 ). 



3.2. Constraints on the physics of dust evolution in 
gas-rich disks 

The /^-values needed for disks larger than 5 AU are sig- 
nificantly lower than the valu es constrained for the ISM 
(« 1.6 — 1.8, see |Draine||2003| ) and are instead consistent 
with the values derived lor disks around PMS stars (< 1, 



e.g. |Natta et al.|2007| [Ricci et al.|2010a| >. So far, the only 
plausible hypothesis proposed to explain these values of 
(3 is that most of the thermal emission observed in the 
sub-mm is coming fro m mm- sized grains, or larger, fro m 
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2012). 



isks can 



the disk outer regions ( |Draine||2006} 

State-of-the-art models of dust evolution in c 
calculate the size-dependent evolution of solids account- 
ing for a variety of different mechani sms, e.g. coagula- 



tion, fragmentation, r adial migration ( Brauer et al.| 2008 
|Birnstiel et al.|[2010a[ ). Sub-mm observations providing 
information on the grain size distribution of dust in disks 
around PMS stars h ave been used to test the predictions 



from these models (Birnstiel et al. 2010b Pinilla et al. 



2012). In these models, the efficiency lor i disk to grow 



solids critically depends on its physical conditions. If 
compared with more massive disks around PMS stars, 
these models predict that at the conditions of BD disks 
(e.g. lower densities, higher relative velocities due to 
radial drift) the growth to mm-grains should be much 
less efficient. Furthermore, in BD disks inward radial 
drift of ~mm-grains is expected to be faster than in T 
Tauri disks by a factor of a few. Therefore, the dis- 
covery of mm-sized grains in the outer regions of BD 
disks, as suggested in the case of p-Oph 102, severely 
challenges models of early evolution of solids. Possible 
ways to explain the presence of mm-sized grains in the 
outer regions of young BD disks are discussed in Pinilla 
et al. (2012, in prep.). These include very low levels of 
turbulence in the disk to decrease the collisional speed 
between solids, and require the action of some physical 
mechanisms which can efficiently halt the inward radial 
migration of these particles, e.g. gravitati onal clumping 
favor ed by the magnetic field in the disk ( Johansen et al. 
[2007| ). 

4. SUMMARY 

We presented ALMA observations of the young Brown 
Dwarf p-Oph 102 at about 0.89 mm and 3.2 mm. We re- 
ported the detection of the disk surrounding p-Oph 102 
in both the dust continuum bands and in CO(J = 3 — 2). 
Dust continuum emission is spatially unresolved, and this 
provides an upper limit on the radius of the dusty disk 
of ~ 40 AU. We derived a disk mass of ~ 0.3 - 1% 
of the mass of the Brown Dwarf, similar to the typi- 
cal values found for disks around PMS stars. The de- 
tection of CO molecular gas emission at the location of 
the Brown Dwarf indicates a gas-rich disk as typically 
found for disks surrounding more massive young stars. 
We measured a spectral index of the continuum emission 
ao. 89-3. 2mm = 2.29 ± 0.16, which is significantly lower 
than the value found in the ISM. By modeling the dust 
emission and from the properties of the CO molecular gas 
emission we constrained the variation of the dust opacity 
with frequency, which is a diagnostic of grain growth in 
the disk. Our result indicates that mm-grains are found 
not only in the outer regions of T Tauri disks, but even 
in the low density environments of Brown Dwarf disks, 
and this challenges our current understanding of dust 
evolution in gas-rich disks. 
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Figure 1. Continuum maps of p-Oph 102. Left panel) Continuum map at 0.89 mm. White contour lines are drawn at —3,3,6,9, 18a, 
where a — 0.22 mJy/beam is the rms noise measured on the map. Right) Continuum map at 3.2 mm. White contour lines are drawn at 
3 and 6cr, where a — 0.031 mJy/beam. In each panel, the yellow filled ellipse in the lower left corner indicates the size of the synthesized 
beam, i.e. FWHM = 0.71" x 0.54", PA = 100 deg at 0.89 mm, and FWHM = 1.82" x 1.50", PA = 66 deg at 3.2 mm. For both maps, a 
Briggs weighting with robust parameter = 2 (natural weighting) was used to maximize the signal-to-noise ratio. 




Figure 2. CO(J = 3 — 2) maps of p-Oph 102. Left panel) Moment CO(J = 3 — 2) map. White contour lines are drawn at 2, 4, 6cr, 
where a = 45 mJy/beam-km/s is the measured rms-noise. Right) Moment 1 CO(J = 3 — 2) map. In each panel, the magenta filled ellipse 
in the lower left corner indicates the size of the synthesized beam, i.e. FWHM = 0.57" x 0.46", PA = 99 deg. For the imaging, we adopted 
a Briggs weighting with robust parameter = and considered only projected baselines longer than 70 kA to highlight the emission from 
structures with angular scales < 3". White pixels have been masked out for the computation of the moment maps. 
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Figure 3. Flux at 1 mm vs spectral index between 1 and 3 mm for disks around single PMS stars and Brown Dwarfs. Differe nt colors 
and symbo ls refer to different stellar /sub-st ellar spectra l types a nd regions as indicated in the plot. Data for Taurus disks are from |Ricci et| 
al. ( 2010a] |2012] ), for Ophiuchus disks from [Ricci et al.| (|2010b| ), for 2M0444 from |Bouy et al | ( |2008| > and Ricci et al. (2012, in prep. J, and 
tor p-Uph 1U2 rrom this Letter. Note that ior same disks the values of the 1 mm riux density has been derived by interpolating between 
nearby wavelengths. The typical uncertainties of the data are shown in the lower left corner. 
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Figure 4. Sub-millimeter flux density vs spectral index for disk models for p-Oph 102. The black dot shows our ALMA data for p-Oph 102. 
Each line represents the prediction of disk models with the same disk outer radius and dust opacity spectral index (3, but increasing disk 
mass from left to right. The values considered for the disk radii and (3 are indicated. Each model was computed by assuming a radial profile 
of the surface density X oc r~ p , with p = 1, and a disk inclination of 70 degrees. For th e sub- stellar prope rties of p-Oph 102, we took a mass 
M* = 0.06 Mq, luminosity L* = 0.08 Lq, effective temperature T e ^ = 2700 K (from Natta et al. 20021). For the computation of the dust 
opacities we considered porous spheric al grains made of astronomical silicate s, carbonaceous materials and water ice (optical constants for 
individual components from Weingartner & Draine 2001 Zubko et al. 1996 Warren 1984 respectively) and adopted a simplified version 
of the fractional abundances used by iFollack et al. (1994), as done m |Kiccret al.| (2UlUa b). Grain sizes are distributed as a power-law 
n(a) oc a~ q with a = 3.0 and betw een a minimum grain size of 0.1 /^m and a maximum grain size which determines the value of (3 (for 
more details, see Ricc i et al.|2010a| ). 



